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ABSTRACT 

Context. The Neel IRAM KIDs Array (NIKA) is a fully-integrated measurement system based on kinetic inductance detectors (KIDs) 
currently being developed for millimeter wave astronomy. In a first technical lun, NIKA was successfully tested in 2009 at the 
Institute for Millimetric Radio Astronomy (IRAM) 30-meter telescope at Pico Veleta, Spain. This prototype consisted of a 27-42 
pixel camera imaging at 150 GHz. Subsequently, an improved system has been developed and tested in October 2010 at the Pico 
Veleta telescope. The instrument upgrades included dual-band optics allowing simultaneous imaging at 150 GHz and 220 GHz, faster 
sampling electronics enabling synchronous measurement of up to 112 pixels per measurement band, improved single-pixel sensitivity, 
and the fabrication of a sky simulator to replicate conditions present at the telescope. 

Aims. The primary objectives of this campaign were: to demonstrate imaging with the new dual-band optics, to evaluate the feasibility 
of hundred-pixel arrays of KIDs, to validate the sky simulator's ability to reproduce real observing conditions, to identify performance- 
limiting noise sources, and to image calibration and astronomically-relevant sources. 

Methods. The imaging sensors consisted of two spatially-separated arrays of KIDs. The first array, mounted on the 150 GHz branch, 
was composed of 144 lumped-element KIDs. The second aiTay, mounted on the 220 GHz branch, consisted of 256 antenna-coupled 
KIDs with a distributed resonator geometry. Each of the arrays was sensitive to a single polarization; the band splitting was achieved 
by using a grid polarizer. The optics and sensors were mounted in a custom dilution cryostat, with an operating temperature of ~70 
mK. Electronic read-out was realized using frequency multiplexing and a transmission line geometry consisting of a coaxial cable 
connected in series with the sensor array and a low-noise 4 K amplifier. 

Results. The new dual-band NIKA was successfully tested in October 2010, performing in-line with sky simulator predictions. 
Initially the sources targeted during the 2009 run were re-imaged, verifying the improved system performance. An optical NEP was 
then calculated to be around 2 • 10 W/Hz''-. This improvement in comparison with the 2009 run verifies that NIKA is approaching 
the target sensitivity for photon-noise limited ground-based detectors. Taking advantage of the larger arrays and increased sensitivity, 
a number of scientifically-relevant faint and extended objects were then imaged including the Galactic Center SgrB2(FIRl), the radio 
galaxy Cygnus A and the NGC1068 Seyfert galaxy. These targets were all observed simultaneously in the 150 GHz and 220 GHz 
atmospheric windows. 
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1. Introduction 

The importance of millimeter and submillimeter astronomy is 
now well established. In particular, three main areas of millime- 
ter continuum research have motivated the rapid development of 
new technologies: 

1. The study of star forming regions in the Galaxy 
(Ward-Thompson 20071). The pre-stellar phases in molecu- 
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lar clouds are hidden by cold dust (around 10 K) which can 
only be observed at submillimeter wavelengths. 

2. The investigation of high redshift galaxies 
(|Lagache et al. 2005 1. The redshift effect at submil- 
limeter wavelengths counteracts the distance dimming 
(Blain et al. 2002l i. 

3. The measurement of cosmic microwave background (CMB) 
temperature anisotropics (either primordial or secondary). 
At a temperature of 2.725 K, the CMB spectrum peaks 
at millimeter wavelengths. Of particular interest is the 
Sunyaev-Zel'dovich effect dBirkinshaw 19991 1. which dis- 
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torts the CMB spectrum at millimeter and radio wavelengths 
and can be used to map the distribution of hot gas in clusters 
of galaxies. 

Throughout the previous decade, instruments utilizing hun- 
dreds of individual bolometers in focal plane arrays have dom- 
inated continuum submillimeter and millimeter astronomy (e.g. 
MAMB02, BOLOCAM). Full-sampling arrays with up to thou- 
sands of pixels in a single array are now reaching maturity, offer- 
ing increased mapping speed and decreased per-pixel manufac- 
turing costs (e.g. Apex-SZ, SPT, SCUBA2, LABOCA). Despite 
these considerable advances, further array scaling is strongly 
limited by the multiplexing factor of the readout electronics. 

A promising alternative to traditional bolometers is the ki- 
netic inductance detector (KID). First demonstrated less than 10 
years ago ( [Day et al. 2003) 1, a KID consists of a high-quality su- 
perconducting resonant circuit electromagnetic ally coupled to 
a transmission line. Typically, a KID is designed to resonate 
from 1 to 10 GHz and exhibits a loaded quality factor exceed- 
ing Ql > 10^. Thus each KID occupies a bandwidth of order 
A/ = //fii -10-100 kHz. A single KID only loads the trans- 
mission line within A/ around its resonant frequency. The KID 
resonant frequency can easily be controlled geometrically during 
the circuit design. It is therefore possible to couple a large num- 
ber of KIDs to a single transmission line without interference as 
long as the inter-resonance frequency spacing exceeds 2A/. The 
inter-resonance frequency spacing and the total bandwidth of the 
measurement electronics sets a limit on the total number of pix- 
els which can be simultaneously read out on a single measure- 
ment line dMazin 2004 ). Current frequency-multiplexing mea- 
surement electronics are limited to at most a few hundred pix- 
els per measurement cable. This is expected to grow into the 
thousands of pixels per cable within the next decade, vastly im- 
proving the multiplexing factor in comparison with existing ar- 
ray technologies. 

In order to absorb incident radiation, it is necessary to 
impedance match the KID to free space at the target wave- 
length. Currently, there are two methods for meeting this cri- 
teria. The first is to use a geometry known as the Lumped- 
Element KID (LEKID) which separates the resonator into an 
inductive meander section and a capacitor At the target wave- 
length, the inductive meander approximates a solid absorber By 
accounting for the resonator and substrate impedances and the 
cavity formed with the sample holder, it is possible to directly 
impedance match the LEKID to free space ( [Doyle et al. 2009[ ). 
Alternatively, a lens and antenna structure can be used to adapt 
the resonator to free space. Due to the increased geometric com- 
plexity, it has proven more difficult to achieve satisfactory device 
performance with antenna structures than with direct absorbtion 
LEKIDs. Despite this drawback, antenna-coupled KIDs remain 
a very active area of research due to their frequency selectivity 
dSchlaerth et al. 2008l l. 

For both LEKIDs and antenna-coupled KIDs, detection is 
achieved in the same manner once the incident radiation has been 
absorbed. In a superconductor, the conduction electrons are con- 
densed into charge carrying Cooper pairs. Mediated by lattice 
vibrations, superconductivity results in an energy gap in the car- 
rier density of states. Incident photons with an energy exceeding 
the gap energy can break a Cooper pair, producing two quasipar- 
ticles and a concurrent change in the complex impedance. The 
result is a shift in the KID resonance frequency which can be 
read out by the measurement electronics. 

We are currently developing a fully-integrated measurement 
system based on KIDs known as the Neel IRAM KIDs Array 



(NIKA). The two primary goals of NIKA are to asses the 
viability of KIDs for terrestrial astronomy and to develop a 
filled-array, dual-band resident instrument for the Institute for 
Millimetric Radio Astronomy (IRAM) 30-meter telescope at 
Pico Veleta, Spain. Based on a custom-designed dilution cryo- 
stat with a base temperature of -70 mK ( IBenoit et al. 20081) . 
a first generation single-band NIKA prototype was previ- 
ously tested at the Pico Veleta telescope in October 2009 
(Monfar dini et al. 2 010). This successful measurement was the 
first to directly compare the performance of LEKID and antenna- 
coupled KID designs. 

Leveraging the experience gained from the first generation 
NIKA, an improved instrument has been designed and tested at 
the IRAM 30-meter telescope in October 2010. This second gen- 
eration system includes a large number of enhancements. Dual- 
band optics, integrating a polarization-sensitive splitter and a 
new baffling structure, allow simultaneous imaging at 150 GHz 
and 220 GHz. Resonator design modifications resulted in im- 
proved single-pixel sensitivity. Faster digital-signal-processing 
electronics enable synchronous measurement of up to 112 pix- 
els for each measurement band. Fabrication of a sky simulator to 
replicate typical measurement conditions at the telescope facil- 
itated improved array testing and quality control. Along with a 
detailed discussion of these system upgrades, we present here 
the results of the October 2010 measurement campaign. This 
includes a discussion of the limiting noise sources, an analysis 
of the system performance determined using calibration sources 
such as planets, and an estimate of the improved full-system sen- 
sitivity. Finally we present astronomically-relevant observations 
of a number of faint and extended sources in both measurement 
bands which were previously unattainable with the less sensitive 
first generation NIKA. 



2. The dual-band optics 




Array 

Fig. 1. NIKA optical design. Two flat mirrors (M6,M7) orient 
the beam on the bi-conic mirror MS. The rays enter the cryo- 
stat through the window, coincident with lens LI. A grid polar- 
izer located at 70 mK between lenses L2 and L3 splits the beam 
and redirects the radiation to two polarization-sensitive arrays, 
the first optimized for 150 GHz and the second for 220 GHz. 
The band-defining filters centered at 150 GHz and 225 GHz are 
mounted dkectly on the flat face of each L3 lens. 



Monfardini et. al.: The new NIKA: A dual-band millimeter- wave kinetic inductance camera for the IRAM 30-meter telescope 



3 



The new dual-band NIKA is engineered to fit the receiver 
cabin of the IRAM telescope in Pico Veleta, Spain. The 30-meter 
primary mirror (Ml) and the hyperbolic secondary (M2, D-2m) 
are installed directly on a large alt-azimuth mounting. The in- 
cident beam is directed into the receiver cabin through a hole 
in Ml using a standard Cassegrain configuration. A rotating ter- 
tiary (M3) provides a fixed focal plane (Nasmyth focus). The 
optical axis, in order to conform to the dimensions of the cabin, 
is deviated by two flat mirrors (M4 and M5). M4 can rotate be- 
tween two fixed positions, selecting either heterodyne or contin- 
uum instruments. 

NIKA re-images the large telescope focal plane onto the 
small sensitive area covered by the KIDs. The demagnification 
factor is around 6.6, achieving a well-adapted scale of 5 arc- 
seconds/mm on the detector plane. This is accomplished us- 
ing two flat mirrors (M6, M7), one bi-conical mirror (M8) and 
three high-density polyethylene (HDPE) corrugated lenses (LI, 
L2, L3). The lens corrugation consists of machined concen- 
tric grooves, providing a soft transition between vacuum (n=l) 
and HDPE (n=1.56) to reduce reflective losses. The size of the 
grooves is 0.4 mm x 0.4 mm in depth and width, while the 
width of the ridges is 0.4 mm. LI is located at room tempera- 
ture and coincides with the cryostat vacuum isolation window. 
L2 is mounted on the screen at 4 K, while the final L3 lenses 
(one per array) are installed at the coldest stage just in front of 
the arrays (~70 mK). A simple grid polarizer with a grid pitch of 
4 fim was inserted at 45 degrees with respect to the main optical 
axis before the final L3 lenses. Since all KID designs currently 
employed in NIKA are sensitive to a single polarization, this is 
an efficient way to realize beam splitting for the two frequency 
bands. A schematic of the optical design is presented in Fig.[T] 

For both arrays, the lenses are telecentric. That is, each point 
of the detector plane is illuminated with the same aperture with 
the chief ray perpendicular to the surface. This strongly reduces 
potential systematic errors related to overlap of the incident 
beam with the intrinsic pixel lobe. The detector plane is over- 
sampled with respect to the diffraction spot size: the full width 
at half maximum is about 3 mm for the 150 GHz band and 2 mm 
for the 220 GHz band. The pixel pitches are 2.25 mm and 1.6 
mm, for final aperture ratios of 0.75 -FA and 0.80 F/1 respectively. 
The useful, projected field-of-view is approximately 1 .6 x 1 .6 ar- 
cminutes for the 150 GHz array and 1x1 arcminutes for the 220 
GHz array. 




350 



Frequency (GHz) 

Fig. 2. Normalized spectral response of the NIKA detectors. The 
3 dB measured bandwidths are 135-169 GHz and 198-238 GHz. 

A series of filters and baffiing are used to reject unwanted ra- 
diation. To reduce off'-axis radiation, a multi-stage, black-coated 
baffle is installed at 4 K between LI and L2. This is further aug- 



mented with a cold pupil at 70 mK. Three infrared-blocking fil- 
ters are mounted on the first two radiative screens held at a tem- 
perature of ~150 K and ~70 K by cold helium vapors. Three ad- 
ditional low-pass filters are installed on the 4 K and 1 K screens. 
The radiation entering the 70 mK stage is thus restricted to fre- 
quencies V <300 GHz. The final band definition is achieved us- 
ing a series combination of a high-pass and a low-pass filter be- 
hind the final lens L3, just in front of the detector array. The 
spectral response of NIKA was characterized using a Martin- 
Puplett interferometer (MPI) (IDurand 20071 1. The results, dis- 
played in Fig. |2] exhibit good agreement with the atmospheric 
transparency windows and coincide with the expected response 
based on the individual filter cutoff frequencies. Taking into con- 
sideration the tabulated HDPE transmission and the individual 
filter specifications, we estimate a total optical transmission co- 
efficient of ~0.4 for the NIKA optics. 

3. The detectors and readout electronics 




Fig. 3. Mounted arrays and pixel micrographs, (a) 144 pixels 
LEKIDs array, (b) 256 pixels antenna-coupled array. KID mi- 
crographs: (c) LEKID and (d) Antenna-coupled KID. For both 
micrographs, light regions indicate where metal is present and 
dark areas where the substrate dielectric is exposed. 



The second generation NIKA implements a 144 pixel ar- 
ray of LEKIDs for 150 GHz detection and a 256 pixel array of 
antenna-coupled KIDs for 220 GHz sensing. An image of the 
mounted arrays is shown, along with micrographs of the indi- 
vidual pixels, in Fig.|3] For the LEKIDs array, fabrication com- 
menced with an argon plasma surface treatment of a 300 fim 
thick high-resistivity silicon wafer (>5 kQ/cm). Next, a 20 nm 
aluminum film was deposited via DC sputtering. The final ar- 
ray structure was subsequently defined with UV lithography fol- 
lowed by wet etching. 

For the antenna-coupled KIDs, a hybrid material structure 
necessitated a more complex fabrication strategy. The basic KID 
consists of a quarter- wave coplanar waveguide (CPW) resonator 
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One end of the resonator is coupled to the measurement trans- 
mission line while the other end is grounded. A twin-slot antenna 
focuses incident radiation on the grounded end. Half the length 
of the resonator center strip, on the grounded end, is made of alu- 
minum; the rest of the center strip on the coupler side, and the en- 
tire ground plane, is made of NbTiN. Fabrication consisted of a 
300 nm thick sputter deposition of NbTiN (IBarends et al. 2010] l 
on a HF-passivated high-resistivity silicon wafer (>10 kQ/cm). 
A 100 nm aluminum film was then sputter deposited through a 
lift-off mask into the plasma-etched slots in the ground plane to 
complete the resonator center strip. Aluminum air-bridges were 
then used to connect the ground planes of the CPW in order to 
suppress microwave mode conversion, reduce the cross talk and 
improve the beam pattern. The radiation is concentrated and fo- 
cused on the antennas by silicon elliptical lens segments of 1 .6 
mm diameter, arranged in a rectangular grid of 1.6 mm spacing. 
This lens array is glued on the back side of the KID sample. The 
misalignment between the lens and the twin slot antenna is less 
than 10 jum, ensuring good coupUng efficiency. 

The general principles of frequency-multiplexed 
KID readout have been described in detail elsewhere 
dSwenson et al. 20091 r^teseFal . 2009 Monfa rdini et al. 20101 1. 
The current digital electronics used for the NIKA readout were 
developed in the context of an international collaboration known 
as the Open Source Readout (OSR) (IDuan et al. 20101 1. This 
system was based on a digital platform known as ROACH, 
itself having been developed in the context of another col- 
laboration known as Center for Astronomy Signal Processing 
and Electronics Research (CASPER) (IParsons et al. 20061 1. 
The ROACH hardware provides powerful signal processing 
capabilities by integrating a field-programmable gate array, an 
on-board power pc, and a variety of high-speed communication 
interfaces. Building on this, the OSR developed new high-speed, 
dual-input analog-to-digital and dual-output digital-to-analog 
interface cards. For NIKA, both of these cards were clocked 
by the same rubidium-referenced external clock generator 
Operating at 466 megasamples per second, the resulting useful 
IF measurement bandwidth of the NIKA readout was 233 MHz. 
In order to drive the individual pixels and subsequently read 
out their state, the NIKA collaboration developed customized 
software to use with the OSR hardware. Similar to standard 
lock-in techniques, the implemented algorithm allowed 112 
separate measurement tones to be generated and simultaneously 
monitored within the IF measurement bandwidth. 

The NIKA readout uses a standard up-down converter con- 
figuration based on two IQ mixers per board to transpose the 
generated IF frequency comb to the resonator operating frequen- 
cies. One ROACH board and a set of IQ mixers was used for 
each array. Currently, the LEKID array operates in the frequency 
range 1.27-1.45 GHz. Within this bandwidth, 104 pixels and 8 
off-resonance blind tones could be used for measurement. The 
antenna-coupled KID array has a central pixel core operating at 
5-5.2 GHz. Due to the larger inter-resonator frequency spacing 
of this array, only 72 core pixels could be simultaneously mea- 
sured out of the total 256 pixels in the array. 

The response of every pixel in an array is measured simul- 
taneously and broadcast via UDP packets by the ROACH elec- 
tronics to the control computers at a rate of 22 Hz. The individ- 
ual pixel responses are composed of a pair of in-phase (/) and 
quadrature (Q) values which result from the final stage of digital 
mixing and low-pass filtering. Theses values can be translated 
into the traditional transmission phase and amplitude using the 
identities 6 = arctan(2//) and amplitude^ - + Q^- An alter- 
native approach is to plot these values in the complex plane. An 
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Fig. 4. KID calibration, (a) I-Q frequency scan across a reso- 
nance. Complex phase is calculated by fitting each resonance 
to a circle (shown); the phase angle is then determined with re- 
spect to the center of the resonance circle as indicated by the 
X. The minimum transmission is marked with a small circle, (b) 
Frequency shift versus complex phase. The x on the curve indi- 
cates the maximum frequency shift when Mars transits the pixel. 



example of a standard frequency sweep around a resonance is 
provided in Fig. Ufa). Small changes in illumination result pri- 
marily in motion around this curve and thus it is convenient to 
define a new angle (p about the center of curvature (7^, Qc)'- 



arc tan 



Q-Qc 
I -I, 



(1) 



where 4>q rotates the plane such that the curve intersects the x 
axis at the resonance frequency /q. 
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Fig. 5. Readout power and frequency optimization. The veloc- 
ity magnitude [(AI/Aff + (AQ/Aff]^'^ (solid, blue) is plot- 
ted along with the optical response [(AI/ATf + [AQIATf]^!^ 
(dashed, red) for each point around a resonance. The velocity 
was measured while making linear frequency steps A/ = 6 kHz 
around the resonance. The optical response was measured by 
sweeping the resonance in frequency with a cold focal plane 
{T ~ 70) K and then at increased temperature AT" -1 K. 

We have observed that the KID response to radiation depends 
critically on the driving power and frequency of the readout elec- 
tronics. While still not fully understood, it is reasonable to con- 
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jecture that this dependence could be caused by: non-linear thin 
film effects, quasiparticle generation or suppression in the su- 
perconductor, or a field-dependent dielectric behavior A plot of 
the velocity around the resonance curve for a linear frequency 
sweep and the measured response to a change in the focal plane 
temperature is shown in Fig. |5]for increasing readout powers. 
For low readout power (< 90 dBm), the responsivity of the KID 
is weak. As the power is increased, the KID responsivity and the 
velocity are seen to increase and to shift to higher frequencies. 
Beyond a critical readout power, a discontinuity in the resonance 
emerges and the frequency shift changes directions and moves 
back toward lower frequencies. While the greatest responsivity 
is obtained in this region, the region near the discontinuity is 
unstable and not currently well understood. For the current mea- 
surement a readout power of ~ -77 dBm, just below the appear- 
ance of the discontinuity, was used. For the frequency selection, 
it is also important to note that while the maxima in responsiv- 
ity and velocity coincide, these generally do not occur near an 
extrema of the traditional 5 21 amplitude or phase or their corre- 
sponding derivatives. The exact relationship between the respon- 
sivity and the shape of the 521 curve is observed to depend on 
a number of implementation details, including pixel cross-talk 
and impedance loading of the transmission line. For this reason, 
the NIKA software now implements a calculation of the velocity 
around the resonance during calibration. 

For low intensity radiation, the density of photo-generated 
quasiparticles in a KID is proportional to the incident photon 
flux. To first order, this results in a linear shift in the kinetic 
inductance and, for thin films, a subsequent linear shift in the 
resonance frequency (ISwenson et al. 20T0l . That is: 

<5/o = -C/^SLk cc -'JLdP, (2) 

where C is a constant, n j is the Cooper pair density, and 6Pi is the 
incident power A plot of the phase (p versus the frequency shift 
from resonance d fo for a typical KID is shown in Fig.|4fb) along 
with the maximum frequency shift during a transit of Mars. 
From this plot it is clear that d(p/df is approximately linear for 
relevant astronomical signals. Thus oc dP, with the constant 
of proportionality being determined during a calibration scan. 
One caveat is that under large background changes, which can 
result from weather changes or telescope repointing, n , changes 
substantially and causes a significant shift in /o. For this reason, 
an automatic procedure was implemented to properly determine 
(Ic, Qc) and 0o before every scan taken on the telescope to as- 
sure that the measurement would remain in the dynamic range 
of the KIDs. This recalibration is currently achieved in under 60 
seconds. 

4. The sky simulator 

In order to replicate real observing conditions and to properly es- 
timate the amount of stray-light on the detectors, we have built 
a testing tool to complement MP interferometer measurements 
and classical chopper tests alternating between hot and cold 
sources. Simply called the sky simulator, the basic idea consists 
in cooling down a large, black disk with the same dimensions 
as the telescope focal plane. This cold disk simulates the back- 
ground temperature in ordinary ground-based observing condi- 
tions. On a telescope the main contributions to the background 
are the atmospheric residual opacity, which is weather depen- 
dent, and the emissivities of the mirrors. The typical background 
temperature of the IRAM Pico Veleta telescope is in the range 



30-100 K. For the sky simulator, the cooUng of the background 
disk is achieved by using a single stage pulse-tube refrigerator 
A large window in the sky simulator cryostat is fabricated in 4 
cm thick HDPE which is sufficiently transparent at millimeter 
wavelengths. The minimum background temperature that can be 
achieved is 50 K, limited by the radiation absorbed by the large, 
black, cold disk. An image of the sky simulator can be seen in 
Fig.|6] 

To simulate an astronomical source, a high-emissivity ball 
with a diameter of 5-10 mm was placed in front of the sky sim- 
ulator window at room temperature. The ball is mounted on a 
motorized XY stage enabling movement with respect to the fixed 
disk. The angular speed on the sky of a typical on-the-fly scan at 
the IRAM telescope is 10 arcseconds/s. Accounting for the -300 
m effective focal length of Ml andM2, this corresponds to a scan 
speed of 15 mm/s on the telescope focal plane. This velocity is 
well within the capabilities of the XY translator A comparison 
of two single-pixel transits, the first taken with the telescope and 
the second with the sky simulator, is shown in Fig. |6jc). 




Time (s) 

Fig. 6. Sky simulator (a) HDPE cryostat window removed. The 
cold disk simulating the sky background is shown; its diameter 
is 24 cm. (b) Fully assembled system. The translation arm, nylon 
wire and high emissivity ball are outlined for visibility. The ball 
moves at a controlled rate in front of the HDPE window, faking 
a real source on the sky. In particular, on-the-fly telescope scans 
are easily simulated jBideaud 20101 1. (c) Sky simulator valida- 
tion. A real on-the-fly scan of Uranus taken at the 30-meter tele- 
scope is shown (red, upper) along with a laboratory on-the-fly 
scan using the sky simulator (blue, lower). The telescope angu- 
lar scan speed was 10 arcseconds/s. The linear scan speed of the 
sky simulator was 30 mm/s, corresponding to a telescope angu- 
lar speed of 20 arcseconds/s. For both traces, the transit of the 
source was observed with a single pixel. 

The sky simulator temperature can be adjusted continuously 
between 50 K and 300 K, allowing an accurate estimation of the 
detector response and a direct determination of the noise equiva- 
lent temperature (NET). To perform this measurement, the spec- 
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tral noise density S„( f), expressed in Hz/Hz ' , is calculated at 
a fixed background temperature. The background temperature is 
then increased by 6T and the shift in the KID resonance fre- 
quency is measured. This yields the KID frequency response R 
expressed in Hz/K. The NET is then given by: 

Sf(f) 

NET(f) = (3) 

Laboratory tests performed on the LEKID array indicated an 
average NET per pixel of 4 mK/Hz'^^ at a standard representa- 
tive frequency of 1 Hz. Accounting for the optical chain trans- 
mission and including a 50% reduction due to the polarizer, this 
corresponds to an optical noise equivalent power (NEP) of ap- 
proximately 2 ■ 10"'^ W/Hz'/2. The best single-pixel LEKIDs 
have an optical NEP, under realistic 4 to 8 pW loading per pixel, 
around 6 ■ 10 W/Hz'^^. These results indicate that no concep- 
tual limitations exist for using KIDs in the next generation large 
arrays of ground-based mm-wave instruments. 

The sky simulator has also been used to estimate the amount 
of undesired stray light incident on the detector arrays. This 
measurement required two steps. First the pixel frequency shifts 
were recorded while the sky simulator was moved in the optical 
axis direction from its customary position at the telescope fo- 
cal plane up to the cryostat window. Next, the sky simulator was 
returned to the telescope focal plane and the background temper- 
ature was then increased until an equivalent frequency shift was 
affected. This resulted in an estimated stray-light temperature 
of ~ 35 K which corresponds to ~ 4 pW of parasitic power per 
pixel at 150 GHz. The unwanted radiation has thus been reduced 
by more than a factor of two compared with the first generation 
NIKA and is now comparable to the best sky conditions at Pico 
Veleta. 

5. Results 

The dual -band NIKA run took place in October 2010. The in- 
strument was installed in the receiver cabin of the IRAM 30- 
meter telescope at Pico Veleta, Spain, and operated remotely 
from the control room. The cool-down of the instrument was also 
performed remotely, taking approximately 1 8 hours to reach the 
operating temperature of 70 mK. 

Astronomical data from the two arrays are reduced off-line 
with dedicated software. The raw data (1,0) are converted to 
complex phase angle using the closest previous KID calibra- 
tion. Then a conversion to an equivalent frequency shift is done 
with the same calibration using the derivative of the frequency 
with the complex phase at the zero phase, as described in Fig.H] 
Data are thus internally converted to frequencies which are as- 
sumed to be linear with the absorbed photon counts, as in equa- 
tion Eq. (|2]l. After opacity correction, and using Mars as the pri- 
mary calibrator, we obtain that the overall median gain is of 14 
mJy/beam/Hz and 9 mJy/beam/Hz for the 1 .4 and 2 mm (220 
GHz and 150 GHz) channels, with a 30% dispersion. The fo- 
cal plane geometry of each array is measured by using scanning 
maps of planets (see Fig. |7]i. 

The fitted focal plane geometry is found by matching the 
pixel position in the array as measured on the wafer to the mea- 
sured position on planets, by optimizing a simple set of pa- 
rameters: a center, a tilt angle and a scaling expressed in arc- 
seconds/mm. Most detectors are within less than 2 arcseconds 
of their expected position. The beam width is also found from 
planet measurements. Typically the FWHM is 12.4 and 16.7 arc- 
seconds for the two arrays (1.4 and 2 mm respectively, see figure 
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Fig. 7. 80% (Blue) and 50% (Grey) contours of the beams as 
measured on Mars on the 2010-Oct-20. Top: 2 mm array (98 
valid KID), Bottom: 1.4 mm array (62 vaHd KID). 



8) with a dispersion of 1 arcsecond. This is close to the diffrac- 
tion limit for the 2 mm array. Pixelisation effects have not been 
removed from this estimate. This might explain why the 1 .4 mm 
beam is larger by 20% than the diffraction size. 

We were able to observe many types of objects. We here 
show only examples of the results that are currently analyzed. 
Usually, the observing scan mode is a simple on-the-fly map- 
ping, composed of constant-elevation sub-scans, with the tele- 
scope sweeping in azimuth back and forth. The detector data 
are processed in the time domain (calibration and filtering) be- 
fore being projected on a common map, using the focal plane 
geometry as described above. The noise is evaluated at the de- 
tector level by histogram Gaussian fitting outside the sources. 
The noise is then propagated to the map level. The final map is 
obtained by coadding the individual detector maps with an in- 
verse square noise weighting scheme. From known sources and 
reproducibility from scan to scan and from detector to detector, 
the photometric accuracy is estimated to be 30%. New methods 
are being investigated to improve this photometric offline pro- 
cessing accuracy. 

We present SrgB2(FIRl) as an example of extended source, 
Cygnus A as an example of a multiple source, and NGC1068 as 
an example of weak source with a new 2 mm flux measurement. 
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Fig. 8. Caption: Beam contours obtained on Mars with levels at 
10 to 90 percent. Azimuth and elevation slices along the center 
are shown on the sides along with the slice of a 2D Gaussian 
fit. The FWHM is 16.7 (resp. 12.4) arcseconds at 2 mm (resp. 
1.4 mm). 



Fig. |9] shows a dual map of the Galactic Center obtained in 
900 seconds of integration time. This complex region reveals at 
least 3 compact sources. The very center has a flux of 76 ± 3, Jy 
and 17.7 + 0.7, Jy at 1.4 and 2 mm, respectively. 

In Fig. (To] we show Cygnus A (3C405), a well-known ra- 
dio source with two prominent radio lobes. The maps were 
obtained with two scans making a total integration time of 
2200 seconds. The central source flux is of 0.76 + 0.09 Jy (resp. 
0.89 + 0.04 Jy) at 1.4 mm (resp. 2 mm). The flux errors are dom- 
inated by photometric noise and not detector noise. The two ra- 
dio lobes (named AB and DE) have a flux (taken at the posi- 
tions given by [ Wright and Birkinshaw (2004)) of 1.63 + 0.09 Jy 
and 1 .65 + 0.08 Jy at 1 .4 mm. At 2 mm, the measured fluxes 
are 2.47 + 0.13 Jy and 2.36 ± 0.10 Jy. The 1.4 mm fluxes 
can be compared with BIMA 1.3 mm interferometer observa- 

The BIMA fluxes of 



tions by Wright and Birkinshaw (2004) 



the nucleus, the AB and DE radio lobes are 0.48, 0.54, and 
0.97 Jy. The spectral dependance {F 

is 



Wright and Birkinshaw (2004) 



oc V"') as measured by 
recovered on the two radio 
lobes with the NIKA camera. Nevertheless, the NIKA fluxes 
are a factor 1.70 higher. We think that this an angular resolu- 
tion effect. For example, using the 1.1 mm flux measured by 
[Eales, Alexander & Duncan (1989)| with a 19 arcsecond beam 
and applying the above spectral dependence, we expect 0.74, 1 . 1 



and 1 .5 Jy at 1 .4 mm for the flux of the nucleus and the two radio 
lobes which is in agreement with what we obtain. 

Fig. [TT] shows a secure detection of NGC1068, a nearby 
galaxy with an active galactic nucleus. The map is obtained 
with 5 scans and a total integration time of 1260 seconds. Its 
flux is 142 + 25mJy at 1.4 mm and 66 + 3mJy at 2 mm. 
This is the central flux measured with a Gaussian of 12 and 19 
arcseconds, respectively. For this map, a sky noise decorrela- 
tion has been used, which is based on a linear regression with 
the detector signals when off'-source. The North East - South 
West extension at 1mm is larger but aligned with interferomet- 
ric IRAM Plateau de Bure Interferometer (PdBI) observations 



(Krips et al. 2006). The PdBI measured flux for the core and the 
jet is 28 + 2mJy at 231 GHz. This is smaller than the NIKA 
flux measurement which is more consistent with the flux of 
170 + 30mJy measured bv lThronson et al. 19871 indicating the 
presence of a diffuse extended component. This component is 
expected (Hildebrand et al. 19771), based on consideration of the 
far infrared spectral energy distribution. It is likely to come from 
heated dust in the circumnuclear region. 

From deep integration on weak sources and using sky noise 
decorrelation, we are able to derive the eff'ective sensitivity of 
the camera, in the present early state of data reduction and sky 
noise subtraction techniques. We obtain a weak-source flux de- 
tectivity of 450 and 37, mJy • i'^^ at 1.4 and 2 mm respectively. 
The 1 .4 mm detectivity is satisfactory for an initial 1 .4 mm KID 
prototype. The 2 mm detectivity shows a major improvement by 
a factor 3 with the best value obtained in the 2009 NIKA run 
(IMonfardini et al. 2 010). This detectivity is almost at the level 
of the state-of-art APEX-SZ TES detectors (ISchwan et al. 20 1 Oi l 
albeit obtained here with a larger telescope. An NET of 6mK-i'^^ 
is deduced from this value. A sensitivity to the SZ effect can be 
obtained as 3 ■ IQ^^hr^^^ in the y Compton parameter for one 



8 Monfardini et. al.: The new NIKA: A dual-band millimeter- wave kinetic inductance camera for the IRAM 30-meter telescope 



Average Smm Radec Map: CygA 



Jy/be 



Average 2mm Radec Map: NGC1068 



■Jy/be am 





Average 1mm Radec Map: CygA 



1^ 



Jy/beam 



Average 1mm Radec Map: NGC1068 



Jy/be 



£0^ 



(l.r) 




Fig. 10. Maps at 2 and 1 .4 mm of the radio source Cygnus A Fig. 11. NIKA maps of NGC1068 galaxy at 2 and 1 .4 mm wave- 

and its two radio lobes. All the three sources are detected at lengths. The pixel size is 2 arcseconds. The maps have been 

both wavelengths. Contours are at levels of 0.088 Jy/beam (resp. smoothed with an 8 arcsecond Gaussian. Contours are at the lev- 

0.332 Jy/beam) and multiples at 1.4 (resp. 2 mm). A Gaussian els of 9.2 mJy/beam x (2, 3, ...7) for the 2 mm channel and 21, 



smoothing of 8 arcsecond has been applied to the maps. 

beam, although this may be complicated by the SZ extension. 
SZ measurements will be reported in a later publication. 

The individual detector performances were determined using 
noise spectra taken while on-telescope following the same strat- 
egy implemented with the sky simulator A typical noise spec- 
trum, taken with the LEKID array during a standard on-the-fly 
scan, is shown in Fig. [12] The data is fit using a sum of two 
power laws: 

Nf{f)^Af + Bf (4) 

The low-frequency fluctuations seen while observing are due 
to the expected source and sky noise (a - -1.35) while the 
higher frequency fluctuations can be attributed to the detector 
noise (fi - -0.15). The detectors noise is relatively flat, in con- 
trast with previous measurements where a Ng j oc /~"°'^ depen - 
dence was observed (Monfardi ni et al. 2010! IGao et al. 20071 1. 
This seems to indicate a substantial reduction, at least in the case 
of the LEKID array, of the intrinsic frequency noise due to ran- 
dom variations of the effective dielectric constant. For the 220 
GHz (1.4 mm) antenna-coupled KID array, the spectral slope 
does not contradict the cited model and exhibits the expected 
/""^ dependence. We recognize however that dark measure- 
ments would be more suited to draw definitive conclusions con- 
cerning fundamental noise properties of the resonators. 



42, 63 mJy/beam for the 1.4 mm channel. 



Using the frequency response of Mars, a primary calibra- 
tor, and the noise spectra, we calculate a raw detector sensitiv- 
ity of the order of 25 mJy/Hz'^^ for the LEKID aiTay and 180 
mJy/Hz'^^ for the antenna-coupled KID array. Both numbers re- 
fer to a single typical pixel and at 1 Hz. In terms of the optical 
NEP, with reasonable assumptions on the optical transmission 
chain and the pixel geometry, we estimate an average 2.3 ■ lO""" 
W/Hzi/2 for the LEKID array and about 1.5 ■ 10"'^ W/Hz'^^ 
for the antenna-coupled array, both at 1 Hz. These values in- 
clude a 50% factor which accounts for the impinging power re- 
duction due to the polarizer No further corrections have been 
applied to account for the detectors intrinsic optical efficiency. 
Millimeter-wave VNA measurements performed at room tem- 
perature on the LEKID array indicate a good radiation cou- 
pling, with a peak absorption (central frequency) exceeding 80% 
(IRoesch et al. 2010) . Three dimensional EM simulations sug- 
gest about 65% for the antenna-coupled array optical absorption. 
The estimated NEP are, for both arrays, in good agreement with 
the preliminary values obtained using the sky simulator 
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Fig. 12. Frequency noise spectrum taken during a typical on-the- 
fly scan over the G34.3 galactic source. A fit of the function 
in equation (4) yields the contribution to the noise due to low- 
frequency fluctuations from the source and sky (a - -1.35) and 
the contribution at higher frequencies from the KIDs and the 
electronics (/3 - -0.15). 

6. Conclusions and perspectives 

Work on the NIKA instrument commenced in November 2008. 
At that time, no mm-wave LEKID had been designed and fab- 
ricated while substantial research on mm-wave antenna-coupled 
KID designs had been conducted dSchlaerth et al. 20081) . In only 
two years, NIKA has developed rapidly and the 150 GHz LEKID 
array sensitivity is now comparable to the existing IRAM instru- 
ment MAMB02 at Pico Veleta. While the sensitivity of the 220 
GHz antenna-coupled array did not attain the target performance 
during this run, no significant barrier exists toward a substantial 
improvement in the near future. 

For both arrays, further investigation is needed to achieve 
optimal photon-noise limited performance. Currently, a number 
of avenues are being pursued in order to achieve this goal: 

1. New single-pixel geometries are being designed which 
promise to reduce pixel cross-talk, minimize detection vol- 
ume, and absorb both polarizations. 

2. Promising materials such as TiN are being characterized 
which may drastically improve the quantum efficiency of the 
detectors (Led uc et al. 2010l l. 

3. The readout procedure continues to be optimized. A poor 
photometric precision, of the order of 30%, was obtained 
during this run. Development of new models of the under- 
lying superconducting physics is necessary to better under- 
stand the pixel response to incident radiation. A new off-line 
procedure is under development to include the d(p/df non- 
linearity and thus improve the photometric accuracy. Also, 
leveraging experience acquired during previous bolome- 
ter array development, we are implementing a modulation- 
based readout procedure that should simplify and improve 
the photometry by providing the frequency response of ev- 
ery pixel in real-time. 

4. New readout electronics are being developed to further in- 
crease the multiplexing (MUX) factor of each readout ca- 
ble. The current instrument allows 112 pixels to be read out 
simultaneously in 233 MHz of bandwidth. Compared with 
the first version of NIKA, this is a 6 fold improvement in 
the mux factor. While a significant increase, digital electron- 
ics are readily available that can further improve the mux 
factor by more than an order of magnitude. In order to har- 
ness this technology, we are developing a NIKA-dedicated 
readout system called NIKEL (NIKa ELectronics). A first 



generation NIKEL prototype, unused in the current measure- 
ment, achieved 128 channels over a bandwidth of 125 MHz 
(fBourrion et al. 201 1 1. The next generation NIKEL system is 
under development and will achieve more than 256 channels 
over ~400-500 MHz of bandwidth, using 12-bit analog-to- 
digital converters. 

5. While greater magnetic shielding was utilized for the second 
generation NIKA system, magnetic field effects still limited 
the system performance. In particular, the Earth's magnetic 
field was evident when changing the telescope azimuth an- 
gle. Changes in elevations did not strongly affect the cur- 
rent detectors, but potentially limit the performance of higher 
sensitivity instruments. Using a large current-carrying coil of 
wire, it is now possible to mimic the observed effects due to 
the Earth's field in the laboratory. This new testing capabil- 
ity will allow a much more robust magnetic screening to be 
implemented and tested before installation at the telescope. 

6. Improved baffling and filters are being designed to continue 
reducing spurious radiation. Compared with the previous 
NIKA run, the current instrument reduced the stray light by 
more than a factor of two but stray light continues to degrade 
the detector performance. 

7. NIKA will move to a cryogen free system. While not affect- 
ing the system sensitivity, the use of a standard dilution cryo- 
stat is a practical limitation prohibiting long-term installation 
of NIKA at the Pico Veleta telescope. The new cryogen-free 
cryostat has already been fabricated and is currently under- 
going testing. 

The current version of NIKA demonstrates the potential for 
KIDs to operate in large, ground-based mm-wave instruments. 
The NIKA project goal is to be a 6 arcminutes field-of-view, 
dual-band resident instrument at the 30-meter Pico Veleta tele- 
scope. In order to preserve the telescope intrinsic angular resolu- 
tion, a Nyquist sampling of the focal plane is targeted (0.50F/1). 
This requires ~1500 pixels at 150 GHz and ~3000 pixels at 220 
GHz. As an intermediate step to achieving this specification, 
NIKA will be upgraded to a cryogen-free, dual-band instrument 
covering a field-of-view of about 3 arcminutes. Assuming the 
present ai0.75F/l sampling, a 224 pixel array working at 150 
GHz and a 489 pixel array at 220 GHz are required. Assuming 
a MUX factor of 256, only three cold amplifiers and six coaxial 
cables are necessary. A forthcoming measurement, implement- 
ing these system improvements is planned for the near future . 

Acknowledgements. We would like to thank Santiago Navan'o, Juan Penalves, 
Frederic Damour, Caisten Kramer, David John, Juan Luis Santaren, Denise 
Riquelme, Salvador Sanchez, Hans Ungerechts, Robert Zylka and all the IRAM 
staff for the excellent technical support during the run. We also acknowledge the 
technical staff at Institut Neel that have built the NIKA cryostat and paiticipated 
in the electronics development, in particular Henri Rodenas, Gregory Garde, 
Anne Gerardin, Julien Minet and in general the Cryogenics and Electronics 
Groups. This work was supported in part by grant ANR-09-JCJC-0021-01 of the 
French National Research Agency, the Nanosciences Foundation of Grenoble 
and Region Rhone-Alpes (program CIBLE 2009). Part of the travel funds for 
the ran have been provided by the French Ministere des Affaires etrangeres 
et europeennes (PHC Alliance 2010). This work is supported, in the UK, by 
STFC. This research, and in particular A. Baryshev, was supported ERC starting 
Researcher Grant ERC-2009-StG 240602 TFPA. Akira Endo is financially sup- 
ported by NWO (VENI grant 639.041 .023) and the Netherlands Research School 
for Astronomy (NOVA). We acknowledge the cracial contributions of Ben Mazin 
(UCSB), Brano Serfass (Berkeley) and the OSR (Open Source Readout for 
MKIDs) collaboration to the NIKA 2010 electronics. 

References 

Bai'ends, R., Vercrayssen, N., Endo, A., et al. 2010, Applied Physics Letters, 97, 
033507 



10 



Monfardini et. al.: The new NIKA: A dual-band millimeter-wave kinetic inductance camera for the IRAM 30-meter telescope 



Benoit, A., Bideaud, A., Camus. P.. et al. 2008, in Millimeter and Submillimeter 
Detectors and Instrumentation for Astronomy IV, Proc. SPIE, 7020, 702009 

Bideaud, A. 2010, PhD Thesis (in French), Universite Joseph Fourier, Grenoble, 
France 

Birkinshaw, M. 1999, Physics Reports, 310, Issue 2-3, 97 

Blain, A. W., Small, I., Ivison, R.-J., Kneib, J.-R, Frayer, D. T. 2002, Physics 

Reports, 369, Issue 2,111 
Bourrion, O., Bideaud, A., Benoit, A., et al. 2011,'arXiv:1102.1314 
Parsons, A., Backer, D., Chang, C, et al. 2006 in the Proceedings of the Asilomai' 

Conference on Signals and Systems, 203 1 
Day, P K., LeDuc, H. G., Mazin, B. A., Vayonakis, A., & Zmuidzinas, J. 2003, 

Nature, 425,817 

Doyle, S., Mauskopf, P., Zhang, J., et al. 2009, in The Thirteenth International 
Workshop On Low Temperature Detectors, AIP Proc, 1185, 156 

Duan, R., McHugh, S., Serfass, P., et al. 2010, in Millimeter and SubmilHmeter 
Detectors and Instrumentation for Astronomy V, Proc. SPIE, 7741, 7741 IV- 
77411V-10 

Durand, T. 2007, PhD Thesis (in French), Universite Joseph Fourier Grenoble, 
France 

Bales, S. A., Alexander P, and Duncan, W. D. 1989, MNRAS, 240, 817 

Gao, J., Zmuidzinas, J., Mazin, B. A., LeDuc, H. G., & Day, P K. 2007, Applied 

Physics Letters, 90, 102507 
Hildebrand, R. H., Whitcomb. S. E., Winston, R., Stiening, R. R, et al. 1977, The 

Astrophysical Journal, 216, 698 
Krips, M., Eckart, A., Neri, R., et al. 2006. Astronomy & Astrophysics, 446, 

Issue 1, 113 

Lagache, G., Puget, J.-L.; Dole, H. 2005, Annual Review of Astronomy & 

Astrophysics, 43, Issue 1, 727 
Leduc, H. G., Bumble, B., Day, P K. et al. 2010, Applied Physics Letters, 97, 

Issue 10, 102509 

Mazin, B. 2004, PhD Thesis, California Institute of Technology, USA 
Monfardini, A., Swenson, L.J., Bideaud, A., et al. 2010, Astronomy & 

Astrophysics, 521, id.A29 
Roesch, M., Bideaud, A., Benoit, A., et al. 2010, in MiUimeter and 

Submillimeter Detectors and Instrumentation for Astronomy V, Proc. SPIE, 

7741, 77410N-77410N-9 
Schlaerth, J., Vayonakis, A., Day, P., et al. 2008, Journal of Low Temperature 

Physics, 151,684 
Schwan, D., Ade, PA.R., Basu, K., et al. 2010, arXiv: 1008.0342 /I 
Swenson, L. J., Minet, J., Grabovskij, G. J., et al. 2009, in The Thirteenth 

International Workshop On Low Temperature Detectors, AIP Proc, 1185, 

84 

Swenson, L. J., Cmciani, A., Benoit A. et al. 2010, Applied Physics Letters, 96, 
23-3511 

Thronson, H., Walker C. K., Walker, C. E.. Maloney, Ph. 1987, The 

Astrophysical Journal, 318, 645 
Wright, M. C. H., Birkinshaw, M. 2004, The Astrophysical Journal, 614, 115 
Ward-Thompson, D., Andr, Ph., Crutcher, R., et al. 2007, Protostars and Planets 

V (Tucson: University of Aizona Press) 
Yates, S. J. C, Baselmans, J. J. A., Baryshev, A. M., Klein, B., & Glisten, R. 

2009, Applied Physics Letters, 95, 042504 



